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Abstract 

From an analysis of magnetohydrodymic turbulent processes, we conclude that 
the frequent observation of outward propagating Alfvenic fluctuations in the 
solar wind can arise from early stages of in situ turbulent evolution, and 
need not reflect corcnai processes as has hitherto been assumed. 
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Spacecraft obstrvations have indicated that in the interplanetary medium a 
high degree of correlation frequently exists between the fluctuating magnetic 
field S and the fluctuating plasma fluid velocity v. These are "Alfvenlc 
fluctuations" with rms fluctuation levels comparable to the mean magnetic 
field and have nearly aligned v and Finite amplitude Alfvenic fluctua- 

tions (with V s ±S, S measured in units of an Alfven speed) are exact 
solutions to the equations of ideal incompressible magnetohydrodynamics (MHD) 
turbulence’ and are candidates for preferred final states of MnD**. 

Though the evolution of MHD turbulence is complex, it can proceed toward 
final states which can be simply described. It has been recently argued’ that 
MHD turbulence can increase the degree of alignment at large times and that 
this process produces interplanetary Alfvenic fluctuations. Evidence in 
support of "dynamic alignment" has accumulated in the form of closure 
calculations’ and from numerical solutions of the MHD equations.^’* Under 
other conditions, evidence has been presented that MHD turbulence will 

instead evolve toward a minimum-energy state, with a decaying ratio of kinetic 
to magnetic* energy. It is not entirely clear what boundaries separate these 
two regimes and this remains an unresolved issue in MHD turbulence theory.* 

Here we present further evidence for dynamic alignment of v and S in two- 
dimensional MHD turbulence from which we see that Interplanetary fluctuations 
need not represent time asymptotic states. Rather, we argue that inter- 
planetary Alfvenic fluctuations are generated in situ in the early stage of 
dynamic alignment. Our model not only provides an expK^nation for inertial 
range alignment of v and S’ in t'ae solar wind, but also accounts for the fact 
that the sign of this correlation (relative to the direction of the mean 
msignetic field) corresponds to the observed "direction of propagation"; away 
from the sun*. We find support for the viewpoint that MHD fluctuations in the 


solar wind r«setnbl£ turbulence*''^ arising from local stirring of the medium. 

There are at least two important questions involved in the generation of 
Alfvinlc states: First, does dynamic alignment occur and why? Second, if a 

high degree of alignmiiitht initially exists, does it inhibit further development 
of turbulent processes, such as spectral transfer? (Perfectly aligned states 
exhibit no spectral transfer, which has been used to suggest^* that inter- 
planetary fluctuations are not, properly speaking, evolving turbulence.) 

We have investigated these questions numerically using an incompressible 
MHD spectral method code”^'*. The equations solved are, in a familiar 
dimensionless form: 

3v/3t v*7v = -7p + 5*76 + v7*v (1) 

3b/3t + V*75 s S*7V + U7*b. (2) 

p is the mechanical plus magnetic pressure and 7*v = 0 = 7«5. The magnetic 

and mechanical Reynolds numbers are Vm and 1/u, respectively. Periodic 

boundary conditions in the x-y plane are assumed and 3/3z = 0. The Fourier 

coefficients v(I<,t) and 6(I<,t) ere advanced from initial (t=0) values, for 

each wave number by integrating the truncated Fourier transformed versions 

of (1> and (2). Time is scaled to the Alfv4n transit time of unit distance 

corresponding to the initial magnetic field. 

The degree (if alignment is measured by the ratio of two ideal "rugged” 

invariants, the cross helicity, H s <v«S>/2 an(;| the energy E = <v* + b*>/2 

(brackets denote a spatial average). Extremal Vlalues of H./E (±1/2) indicate 

c 

perfect alignment. 

Four simulatio'hs are discussed here, designated A, B, C and D. Runs A, B, 


and C were obtained with a 128^128 code with Reynolds numbers of 400. Includ- 
ed wavevectors have integer components ranging from 1 to 64 and their negative 
counterparts. At t*0, all non-vanishing Fourier amplitudes were confined to a 
circular annulus 9 < <25 and had a fixed kinetic energy, twice as much 

magnetic energy, and phases chosen as follows: The phase of b(i«,0) was chosen 
randomly and a phase angle e was chosen between v(i<,0) and S(i<,0). This 
specifies the modal cross helicity, H (i<) = Re v*(i<) *S(i<) . Run A is 
initially highly aligned with arccos(0.95) s 9 and H./E = 0.45. Run B is only 
moderately aligned with arccos(.3) = 0 and H^/E = 0.14. Run C has the phases 
of 8 and V both randomly chosen so that H a 0. Initially run D (64x64 grid; 
Reynolds numbers 400) consisted of the homogeneous turbulence used in A, B, 
and C superimposed on a large scale shear flow and a correlated large scale b. 
At t=0, H^/E= -0.4. Time steps were either 1/256 or 1/512. Other runs, not 
described here, displayed similar results. 

Direction averaged modal energy spectra for runs A and C are shown in Fig. 
1, both at t=0.19. Initial energy spectra in both runs were identical. 
Spectral transfer out of the initially excited annulus is not noticeably 
inhibited by the high degree of alignment in run A relative to what it is in 
run C. The nearly power-law inertial range in run A occurs at a lower level 
than for the « 0 case C, indicating that initial alignment can at most slow 
transfer, but cannot prevent turbulence from occurring. 

Fig. 2 illustrates the temporal development of dynamic alignment. The 
ratio |H^/E1 for runs A and B increases nearly monotonically, in agreement 
with previously reported results.*’^** 

A physical explanation for dynamic alignment is suggested by rewriting (1) 
and (2) in the Elsasser variables z* = v ± b. We focus on the low-k region, 
neglecting dissipative terms and combine eqs. (1) and (2) into 3z”/3t s 


-z'^'yz^ - 7p. p is determined by using 7*z* = 0 to yield a Poisson equation. 

Consider now the case |z'*’| » |z”i. Note that z* is nearly independent of 
time. f3inc« 7p is of order z”»7z'*’t the fractional rate j 3zV3t|/ |z'*’l goes to 
zero as z” 0. However, the fractional rate l3z”/3tl/|z“l remains of 0(1) as 
z” ♦ 0. The spectrum of z* is effectively frozen, but that of z“ continues 
i.0 evolve and to acquire higher wave number components. In fact, at every 
time step the z” spectrum will spread to all additive combinations of wave 
numbers that comprised those present in z‘*'(0) and z“ at the previous time 
step. The transfer of the z“ spectrum to higher wave numbers will continue at 
a fractional rate that remains finite even as zT becomes small.* Eventually, 
z" will be dissipated as it reaches the higher k values, ^.nhancing the 
inequality |z | » |z t. We believe this mechanism, by wh,*,ch the majority 
species cannibalizes the minority one by sending it to high wave numbers to be 
dissipated, is the essential mechanism involved in dynamic alignment. 

The simulations reflect this qualitative behavior; Fig. 3 shows the 
normalized direction-averaged cross helicity spectrum 2H^(k)/E(k) vs. k at an 
early stage of runs A, B and D. The transfer of excitation to k's higher than 
those initially present consists almost entirely of cross helicity contribu- 
tions of opposite sign td those in the energy-containing low-k part of the 
spectrum. For runs A and B, z“ is transferred to toward the dissipation 
range, as is z* for run D. 

This effect apparently dominates for times of the order of an Alfven 
transit time or an ”eddy turnover time.” Subsequently, the simulations 
indicate that the normalized cross helicity spectra undergo a further 
systematic evolution as illustrated for run A in Fig, 4. The inertial range 
gradually loses its one-signed character as the majority species slowly moves 
out of the low k*s. For highly correlated runs (such as A), minority species 
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domination of the highest computed k*s persists for very long times (beyond 
('^320, not shown). 

In runs with less aligned initial correlation such as B, behavior similar 
to that of in Fig. 4 is observed, bub with less pronounced minority sp<>cies 
domination of high k's (not shown). The high k part of the spectrum can 
become more or less equally populated with and z~. Because and z~ have 
equal dissipation rates, dynamic alignment continues,' leaving an excess of 
majority species in the energy containing scales. In no case have we seen 
broad-band domination of the inertial range by the majority species before 
very long times (> 20). However, the early time appearance of an inertial 
range populated by the minority species and the very slow spread of the 
majority species from energy containing bo higher k are features common bo all 
our simulations. 

These results suggest a mod'l for the turbulent production of interplaneta- 
ry Alfvenic fluctuations. Consider a volume of solar wind plasma containing a 
mean magnetic field lying entirely in a single magnetic sector.** The large 
scale velocity field is radially outward, so that the sign of the radial 
component of the mean magnetic field determines the sign of the large scale 
contribution to H s H (large). H.(large) is generally nonzero in the 

Q C Q 

vicinity of high speed streams (where Alfvenic fluctuations are usually 
observed) even in the reference frame moving with the mean solar wind 
velocity. Now assume that inertial range fluctuations are due, at least in 
part, to nonlinear couplings which drain energy in the large scale fields and 
pump the small scale fluctuations. If the energy transfer produces dynamic 
alignment consistent with the results described above, then at "early times" 
the inertial range of interplanetary fluctuations would contain a sign of . 
cross helicity opposite to that present in the large scale energy containing 
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structured; th« cro^is h<»llcity of th« inertial range fluctuations, say 

H^(lnertial) , would have a sign corresponding to outward propagation. 

This nodal is consistent with spectral analysis of highly Alfvenic periods 
which often shows that long wavelength fluctuations have H.(k) opposite in 
sign to the inertial range fluctuations.^^ These "inward propagating" long 
wavelength waves would be the transition in the spectrum between the large 
scale energy containing structures which are driving the turbulence and the ■ 
inertial range fluctuations which are being driven. 

Durlhg the later stage, as the turbulence gets "older", the inertial range 
may become equally populated with both z* and z" oa envisioned in Ref. 5. One 
might anticipate that interplanetary plasma which has spent a longer time 
being "stirred" is less likely to be Alfvenic. One reported observation at 5 
AU behind a high speed stream is possible evidence of this. In* that spectrum, 
the cross helicity in the inertial range oscillated in sign and was generally 

i 

not close to its maximal values (Ref. 11, Fig. 11), indicating weak inward and 
outward propagating Alfvenic fluctuations. 

Another possibility,’ Is that dynamic alignment enhances a seed correla- 
tion of V and 6 so that by 1 AU highly Alfvenic fluctuations would be seen. 

In contrast to our model, the majority z-field would represent the Alfvenic 
fluctuations. One difficulty with this model is the limited time available 
for turbulent evolution as the solar wind flows between the sun and 1 AU. A 
substantial fractional increase of the v-S correlation requires tens of Alfven 
transit times (cf. Fig. 2). For nominal Interplanetary conditions at 1 AU, 
the energy containing scale is approximately 10’* cm and the Alfven 

speed is approximately one tenth the flow speed of about 4 x 10' cm/s. Thus 
only a few Alfvin transit times elapse for each AU of mean plasma motion and 
it is unlikely that there is time for a substantial degree of dynamic 
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alignmtnt to occur by 1 AU. Furtharmorc, an additional mtchanlsm would be 
required to account for the outward traveling sense of the v-b correlation. 

Our model naturally generates outward propagating fluctuations. The early 
stages of MHD turbulence produce an Inertial range with fluctuating velocity- 
magnetic field correlation opposite in sign to the correlation of the large 
scale velocity and magnetic fields which d* ive the turbulence. This can 
account for the widespread occurrence of Alfv^nic fluctuations in the solar 
wind. The analytical model presented is essentially unchanged by the addition 
of a mean magnetic field and is applicable bo any dimensionality. Simulations 
in two dimensions have shoivn the mechanism to act starting from a variety of 
Initial conditions. 
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Figurt Captions 


Figura 1, Modal total anargy spaotra at t » 0.19 for runs C (uncorrelated 
initial fields) and A (highly correlated Initial fields:) . The Inltal 
spectra were identical except for the values of ({^(k). Spectral transfer 
occurs in both cases. 

Figure 2. The ratio H/E vs. time for runs A and B, showing the near monotone 

C 

increase in the degree of alignment of v and S due to turbulence. 

Figure 3. Normalized cross helicity spectrum for runs A and B at t s 0.19 and 
D at t » 0*4. Data in Figs. 3 and 4 were smoothed by a 9 point moving 
average before plotting. In all three cases, the transfer of cross 
helicity to high k is almost exclusively opposite in sign to that of the 
low k tsO state. 

Figure 4. (a) ^un A at t s a. 98 when transfer of the Initially dominant 

z'^-mode into the Inertial range begins resulting in a H (k) mixed in sign, 
(b) Run A in a later stage (t s 7.8) where the z* field dominates the 
cross helicity spectrum except at the highest k's. 
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